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pq 1 ! ABSTRACT 

Centaurus B is a nearby radio galaxy positioned in the Southern hemisphere close to the Galactic plane. Here we present a detailed 
analysis of about 43 months of accumulated Fermi-LKT data of the y-ray counterpart of the source initially reported in the 2nd Fermi- 
LAT catalog, and of newly acquired Suzaku X-ray data. We confirm its detection at GeV photon energies, and analyze the extension 
and variability of the y-ray source in the LAT dataset, in which it appears as a steady y-ray emitter. The X-ray core of Centaurus B 
is detected as a bright source of a continuum radiation. We do not detect however any diffuse X-ray emission from the known radio 
y~[ \ lobes, with the provided upper limit only marginally consistent with the previously claimed ASCA flux. Two scenarios that connect 

the X-ray and y-ray properties are considered. In the first one, we assume that the diffuse non-thermal X-ray emission component is 
not significantly below the derived Suzaku upper limit. In this case, modeling the inverse-Compton emission shows that the observed 
y-ray flux of the source may in principle be produced within the lobes. This association would imply that efficient in-situ acceleration 
of the radiating electrons is occurring and that the lobes are dominated by the pressure from the relativistic particles. In the second 
scenario, with the diffuse X-ray emission well below the Suzaku upper limits, the lobes in the system are instead dominated by the 
■ magnetic pressure. In this case, the observed y-ray flux is not likely to be produced within the lobes, but instead within the nuclear 

parts of the jet. By means of synchrotron self-Compton modeling we show that this possibility could be consistent with the broad-band 
data collected for the unresolved core of Centaurus B, including the newly derived Suzaku spectrum. 

\& \ Key words, acceleration of particles — radiation mechanisms: non-thermal — X-rays: galaxies — gamma rays: galaxies — galaxies: 

* , jets — galaxies: individual (Centaurus B) 



1 . Introduction UHECRs robustly identifies those AGN as UHECR sources, or 

implies instead that near by AGN simply trace the mass distribu- 
CentaurusB (PKS 1343 -601) is a relatively nearby (z = 0.0129, tion in the Universe (see lMoskalenko et al.ll2009h . 
1 I West & Tareng hi 1989, distance D = 56Mpc) radio-loud active Searches for the members of the galaxy cluster surround- 
H galaxy. Despite its high radio flux — Centaurus B is in fact the ing Centaurus B at optical wavelengths are significantly ham- 
■ fifth brightest radio galaxy in the sky (I Jones et al.ll200ll) — it is pered by a heavy foreground extinction in the direction of the 
relatively poorly studied. This is mainly due to both its proxim- source. Near-infrared surveys showed an enhancement of galax- 
ity to the Galactic plane and the fact that it is not accessible from ies around the discussed object, but not as significant as ex- 
the Northern hemisphere. Initial studies suggested Centaurus B pected, and in addition consisting of mostly late-typ e systems 
was surrounded by a rich cluster of galaxies hidden behind the (Nagavama et al. 2004; Kraan-Korteweg et al. 2005; Schroder 
Milky Way and related to one of the largest concentrations of et al. 120071) . Hence it seems unlikely that Centaurus B resides 
mass in the local Universe, the Great Attractor ( Lvnden-Bell in a particularly massive cluster, although a somewhat less rich 
et al. 1988; Kraan-Korteweg & Woudt 1999; Kraan-Korteweg environment like a poor cluster or a galaxy group remain valid 
& Lahav 2000). In addition, Centaurus B was found close to the options. This conclusion is in agreement with the X-ray stud- 
arrival detection of the ultra-high energy cosmic ray (UHECR) ies presented in the literature, which did not reveal a pres- 
even t detected recently by the P. Auger Observatory (Abrah am ence of any extended atmosphere of hot ionized gas surround- 
et al. l2007l ; lNagar & Matulic h 2008; Moskalenko et al.l l2009). It ing the source. In particular, the Einstein Observatory detected 
is currently not clear whether a positional coincidence between only an X-ray point source coincident with the radio core of 
nearby active galactic nuclei (AGN) and the arrival directions of Centaurus B. Deep observations with the ASCA satellite con- 
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firmed the Einstein detection, and in addition provided some ev- 
idence f or a low surface-bri ghtness diffuse X-ray halo around 
the core (iTashiro et al.|[l99 8). The total luminosity of this halo, 
^2-iokev - 3.7 x 10 41 ergs -1 , which was proposed to be asso- 
ciated with the extended radio lobes in the system, is however 
orders of magnit ude below that exp ected for a rich cluster (see 
the discussi on inlJones et al .1120011) . Using the ROSAT All Sky 
Survey data jEbeling et al.l (120021) did not find any X-ray feature, 
neither point-like nor extended, at the source position, and thus 
excluded large amounts of X-ray emitting gas analogous to those 
found in rich clusters. 

A comprehens ive study of Centau rus B at radio wavelengths 
was presented by iJones et all (120011) . The compiled maps re- 
vealed prominent jets surrounded by diffuse lobes. The largest 
angular size of the radio structure is ^ 24', corresponding to 
a projected size of about ^ 380 kpc for a conversion scale of 
265 pc arcsec -1 . The actual physical extension of the lobes in 
Centaurus B may be larger than this, because of a likely in- 
termediate inclination to the line of sight implied by the jet- 
to-counterjet brightness and polarization asymmetry. Both the 
large-scale radio morphology of the source, as well as its total 
radio power (Li.4ghz - 5 X 10 41 ergs -1 ), suggests a classifica- 
tion of the object as an intermediate case between Fanaroff-Riley 
class I (FR I) and cla ss II (FR II). 

As noted before, Ta shiro et al. I (fl998h claimed the X-ray de- 
tection of the core and the tentative discovery of an extended 
halo coinciding with the radio lobes of Centaurus B. They argued 
that the most likely origin of the diffuse X-ray halo, best fitted by 
a power-law model with the photon index Tx - 1.88 + 0.19, is 
inverse-Compton emission of the non-thermal electrons within 
the lobes interacting with the cosmic microwave background 
(CMB) photon field. This detection, in f act only the second of 
its kind following th e case of Fornax A (iFeigelson et al.|[l995l : 
IKaneda et all 1 1995b , was followed by a number of analogous 
studies targeting other radio galaxie s with the XMM-Newton , 
Suzaku and Chandra X-ray satellites feataoka & Stawarj|2005l : 
ICroston et al.l 120051 ; llsobe et all I2011L and references therein) . 
For Centaurus B, the Chandra observations resolved the inner 
parts of the main (SW) jet at arcsec sc ale thanks to its excellent 
spatial resolution (iMarshall et al.l2005l) . The detection of the SW 
jet on slightly larger scales was also re ported in the prelimi nary 
analysis of the XMM data presented in lTashiro et al.l (120051) . 

Here we present a detailed analysis of Centaurus B at y-ray 
frequencies with the Large Area Telescope (LAT) onboard the 
Fermi satellite, together with an X-ray analysis of the Suzaku 
data we have acquired very recently. Centaurus B is associated 
with the 2FGL J1346.5 -6027 object in the Ferm i -LAI Second 
Sourc e Catalog (2FGL; lAckermann et al.l 120111 : iNolan et all 
120121) . In our investigation we focus on the extended lobes in 
this system, since their high integrated radio flux (dominating 
over that of the unresolved core by about an order of magnitude) 
and exceptional extension on the sky (about half of a degree) 
suggest that Centaurus B could be one of the few nearby radio 
galaxies for which the diffuse lobes may be detected and re- 
solved, or at least viewed as extended, in y-rays. The motivation 
of the study came from the robust identifi cation of several oth er 
radio galaxies with Fermi-L AI sources (lAbdo et al.ll200 9a.b, 
|20 1 Pel : iKataoka et all 120 1 ll) . including the case of Centaurus A 
for which the 600kpc-scale lobes have in fact been suc cessfully 
resolved at GeV photon energies (lAbdo et al.ll2010bllah . These 
studies have provided important clues for understanding energy 
dissipation processes within the extended lobes of radio-loud 
AGN in general, and particle acceleration mechanisms related 
to the production of UHECRs in particular (see in this context 



lHardcastle et al ] l2009al :lQ 7 Sulli van et al ] 120091 : iPe'er & Loebl 
120121) . 



2. Fermi-LKT Data 

The Fermi Gamma-ray Space Telescope was launched in 2008 
June 1 1 . The LAT onboard the Fermi satellite is a pair conversion 
telescope, equipped with solid state silicon trackers and cesium 
iodide calorimeters, sensitive to photons in a very broad energy 
band from ~ 20MeV to > 300 GeV. The LAT has a large ef- 
fective area (~ 8000 cm 2 above 1 GeV for on-axis events), in- 
stantaneously viewing 2.4 sr of the sky with a good angu- 
lar resolution (68% containment radius better than -1° above 
1 GeV). The tracker of the LAT is divided into front and back 
sections. The front section composed of thin converters has a 
better angular resolution (68% containment radius ^ 0.7° at 
1 GeV) than that of the back section (68% containment radius 
^ 1.2° at 1 GeV) composed of thicker converters. Details of 
the LAT instrument and da ta reduction are described in Atwood 
et al. (120091) . 

The LAT data used here were collected during ~ 43 months 
from 2008 August 4 to 2012 February 18. The source event class 
was chosen and photons beyond the Earth zenith angle of 100° 
were excluded to reduce the background from the Earth limb. 
We also applied the ROI-based zenith angle cut which excludes 
times where the region of interest intersects a given zenith angle 
cut. The P7SOURCE_V6 instrument response functions were 
used for the analysis presented in this paper. 

We utilized gtlike in the Science Tools analysis pack- 
ag43 for spectral fits. With gtlike, a binned maximum like- 
lihood fit is performed on the spatial and spectral distribu- 
tions of the observed y-rays to optimize spectral parameters 
of the input model taking into account the energy dependence 
of the point- spread function (PSF). Analysis using gtlike is 
done in a 20° x 20° region around Centaurus B, which is re- 
ferred to as a region of interest (ROI). The fitting model in- 
cludes other point sources whose positions are fixed at the val- 
ues provided in the 2FGL catalog. The Galactic diffuse emission 
is modeled by gal_2yearp7v6_v0.fits while an isotropic com- 
ponent (instrumental and extragalactic diffuse backgrounds) by 
isotropic_p7v6source.txt. Both background models are the stan- 
dard diffuse emission models released by the LAT team. In each 
gtlike run, all point sources within the ROI and diffuse com- 
ponents in the model are fitted with their normalizations being 
free. In the analysis we also include 2FGL sources outside the 
ROI but within 15° of Centaurus B, with their parameters fixed 
at those given in the 2FGL catalog. 

2.1. Detection and Source Localization 

Figure [T] shows a smoothed count map of a 10° x 10° region 
around Centaurus B above 1 GeV, a corresponding background 
model map, and the background- subtracted count map. Here 
only front events are used to achieve the best angular resolution. 
The map includes contributions from the Galactic diffuse emis- 
sion, isotropic diffuse background, and nearby discrete sources. 
The model parameters of diffuse components and nearby sources 
are set at the best-fit values obtained by gtlike using the data 
above 1 GeV. In this fitting model, the 2FGL source associ- 
ated with Centaurus B (2FGLJ1 346.6-6027) is not included. 
The y-ray source associated with Centaurus B is visible in the 



1 available at the Fermi Science Support Center, 
http ://f ermi .gsfc.nasa.gov/ssc 



2 



J. Katsuta et al.: Fermi-L AT and Suzaku Observations of the Radio Galaxy Centaurus B 



2 



o 
O 



[counts / pixel] 
2.5 5.0 7.5 10.0 



[counts/pixel] 
2.5 5.0 7.5 10.0 



[counts/pixel] 
0.00 0.25 0.50 0.75 1.00 



^ — 1 


^ — i 


^ 1 


1 1 1 _jj 1 

+ + 

+ 

+ 

Centaurus B 

" + + + + 

+ +_ 

i i i i i 


i i i _p i 
+ + 

+ 

+ 

Centaurus B 

" + + + + 

l l 1 l l 


1 1 1 _|j 1 
+ + 

+ 

+ 

Centaurus B 

" + + + H- ' 

+ +_ 

i i i i i 



314 312 310 308 306 314 312 310 308 306 314 312 310 308 306 
Galactic longitude Galactic longitude Galactic longitude 



Fig. 1. (Left) Fermi-hKT count map using front events above 1 GeV around Centaurus B in units of counts per pixel. The pixel size is 0!05. 
Smoothing with a Gaussian kernel of cr = (TA5 is applied. The background 2FGL sources contained in the region are shown as green crosses. The 
position of Centaurus B is indicated as a white cross. (Middle) Background model map. (Right) Background- subtracted count map. 




background- subtracted map. The figure also shows residuals on 
the Galactic plane. These may be attributed to the imperfection 
of the Galactic diffuse model and/or the contributions from dis- 
crete sources not resolved from background, since the residuals 
are seen exclusively on the Galactic plane where the Galactic 
diffuse emission components have high intensities and complex 
spatial structures. This interpretation is also supported by the 
fact that all bright residuals are diffuse and/or accompanied by 
nearby 2FGL sources. On the other hand, Figure Q] indicates that 
the source associated with Centaurus B is a discrete feature in the 
residual map. Actually, the region of Centaurus B is relatively off 
the plane: (/, b) = (309°7, 1°7), where the Galactic diffuse emis- 
sion is almost uniform with Galactic longitude and decreases 
monotonically with Galactic latitude. 

To study this y-ray source in detail, we generated a Test 
Statistic (TS) map around Centaurus B. TS is defined as 
-2A ln(likelihood) obtained by gtlike between models of the 
null hypothesis and an alternative. In this paper, we refer to TS 
as a comparison between models without a target source (null 
hypothesis) and with the source (alternative hypothesis) unless 
otherwise mentioned. In the presented TS map, TS at each grid 
position is calculated by placing a point source with a power-law 
energy distribution with photon index set free. 

Figure [2] shows the close-up TS map of Centaurus B, using 
front data above 1 GeV. The figure clearly shows that significant 
discrete y-ray emission is associated with Centaurus B. In addi- 
tion, the figure seems to suggest that the gamma-ray emission 
is elongated in the direction of the radio lobes. To test this indi- 
cation quantitatively, we compare the likelihood of the spectral 
fit using front and back data in the range 0.2 - 200 GeV for a 
point source and an elliptical shape compatible with the extent 
of the radio lobes (22' x 10' in size; see Figure [2]) with a uni- 
form surface brightness. Here similarly a power-law function is 
assumed to model the source spectrum with its normalization 
and index set free during the fitting procedure (see below). The 
improvement of TS in the case of the ellipse (ellipsel) with re- 
spect to the point-source case, TS ex t = TS e i p i - TS psc , is 2.4. 
Figure [2] also indicates an anisotropy of the gamma-ray emis- 
sion. To test this, we rotate ellipsel in 10° steps with its cen- 
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Fig. 2. Fermi-LAI TS map obtained with a maximum likelihood analy- 
sis using front events above 1 GeV in the vicinity of Ce ntaurus B. Linea r 
green contours of the 843 MHz radio map are overlaid (lMcAdamil991h . 
The blue ellipse compatible with the radio lobes is also shown. Black 
cross and contours represent the estimated position and the positional 
errors at 68% and 95% confidence levels in this analysis, respectively. 
Magenta cross and ellipses represent the position and the errors at 68% 
and 95% confidence levels of 2FGL Jl 346.6-6027. The position of the 
core of Centaurus B is represented as a diamond. 



ter fixed and fit the data for each case. The TS in the case of 
ellipsel is the highest, and TS gradually drops with the rota- 
tion to a minimum in the case of the ellipse (ellipsel) which is 
rotated by 90° from ellipsel. The difference between the corre- 
sponding values of TS is TS rot = TS e i p i - TS e i P 2 = 1.4. To test 
the significance of these improvements, we simulate 1000 front 
and back datasets in 0.2 - 200 GeV for the same time period as 
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that of the observation by gtobssim in the Science Tools analy- 
sis package 1 , using a model where a point source for Centaurus 
B, the Galactic diffuse emission, and the isotropic backgrounds 
are included. The point source is set at the best-fit location and 
has the spectral shape of a power-law function (see below). All 
components in the model are set at the best-fit values in the 
range 0.2 - 200 GeV (see § 12.21) . Each simulated dataset is fit- 
ted with models which include Centaurus B, the Galactic diffuse 
emission, and the isotropic backgrounds. The normalization of 
all components and the index of Centaurus B are set free. We 
compare TS between cases that the spatial shape of Centaurus B 
is the point source, ellipse! , or ellipse!. Note that the position 
of the point source is fixed at the position estimated by using 
> 1 GeV front data for each simulated dataset (i.e., the same es- 
timation as that for the observed data). The results show that the 
probability for obtaining TS ex t > 2.4 and TS rot > 1.4 at once is 
3%. This indicates a weak hint of the spatial elongation along the 
radio lobes in the y-ray band, although we can not significantly 
confirm it. Data that will be accumulated in the coming years of 
the Fermi mission will provide us with better spatial information 
to discuss its extent in more detail. 

We estimate the position of the y-ray source to be (a, 
£)=(206?63, -60?44), following the localization of the highest 
TS (TS max = 23) in the TS map (Figure O. The positional 
errors at 68% and 95% confidence levels are about 0?06 and 
0?11, which are delineated by contours of decrements of TS 
from TSmav by amounts of 2.3 and 6.0. respectively ( Mattox 
et al. 119961) . The estimated position with the 68% error con- 
tour is consistent with the position of 2FGL Jl 346.6-6027 at 
(a, £)=(206?66, -60?47), and contains the core of Centaurus B 
at (or, £)=(206?704, -60?408) dJones et alJliOQlh . TS at the es- 
timated position is 29 using front and back data in the range 
0.2 - 200 GeV, corresponding to a 4.5 cr detection. The TS value 
is lower than in the 2FGL catalog (TS ~ 39 for 0.1 - 100 GeV), 
despite the longer integration time here. As we show in § 12.21 
the source is not significantly variable, although its measured 
flux was greatest in the first half of the time interval analyzed 
here, which overlaps most of the 2-year interval analyzed for 
the 2FGL catalog. In addition, the likelihood analysis for the 
2FGL catalog treated front and back events separately, which 
increased the sensitivity, and TS values, slightly. Hereafter, the 
y-ray source associated with Centaurus B is analyzed as a point 
source with this estimated position. 

2.2. Spectrum and Time Variability 

Figure [3] shows the spectrum of the source associated with 
Centaurus B using front and back events. The spectral energy 
distribution (SED) is obtained by dividing the 0.2 - 20 GeV 
energy band into four energy bins. The indices of the other 
sources are fixed at the values given in the 2FGL catalog. The 
Centaurus B source is fitted with a simple power-law function 
in each energy bin with the photon index fixed at 2.6, the value 
obtained from a broadband fitting over the larger 0.2 - 200 GeV 
energy range (see below). In each fit, all the sources within the 
ROI as well as the diffuse components in the model are fitted 
with their normalizations set free. 

We tested a possible spectral steepening at the highest photon 
energies by performing likelihood-ratio tests between a power- 
law (the null hypothesis) and a a log parabola function (the al- 
ternative hypothesis) for the 0.2 - 200 GeV data. The resulting 
Test Statistic is TSlp = -21n(LpL/LLp) = 0.0, which means 
no significant curvature is present in the spectrum. The param- 
eters obtained with the power-law model are the photon in- 
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Fig. 3. Fermi-LAI spectrum of the Centaurus B source. The error bar 
indicates statistical errors of lcr. The arrow represents the upper limit 
on the flux at 95% confidence level. The blue line represents the best-fit 
power-law function from a binned likelihood fit in the 0.2 - 200 GeV 
range. 



dex Y y = 2.6 ± 0.2 and the integrated 0.2 - 200 GeV flux of 
(1.9 ± 0.5) x 10" 8 photon cm" 2 s" 1 . We note that the uncertain- 
ties of the effective area translate to unc ertainties in the flu x 
of < 10% and < 0.1 for the index (lAckermann et al.ll2012h . 
The best-fit power-law function is represented as a blue line in 
Figure [3] We note that the flux and index obtained in the likeli- 
hood analysis change negligibly if the position of 2FGL J 1346.6- 
6027 or the Centaurus B core is used instead for the Centaurus B 
source. 
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Figure |4] shows the light curve of the Centaurus B source us- 
ing front and back data in the 0.2 - 200 GeV range. We divide 
the whole LAT observation period (about 43 months) into only 
four bins, given the low flux of the source. In each time bin, the 
flux of the source is evaluated by performing gtlike with the 
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Table 1. Log of Suzaku observations 



Name 


Observation ID 


Start time (UT) 


End time (UT) 


Net Exposure 1 (ks) 


Centaurus B 


806017010 


2011 Jul 16 20:15:14 


2011 Jul 18 11:50:13 


93.6 


4U 1344-60 


705058010 


2011 Jan 11 19:01:05 


2011 Jan 14 11:24:11 


93.9 



Notes. (t) After the event screenings. 



Galactic diffuse and isotropic components fixed to the best-fit 
values obtained for the whole time range, and the photon index 
fixed at the value 2.6 obtained from a broadband fitting for the 
range 0.2 - 200 GeV. To test the variability quantitatively, we 
calculate each value of the log likelihood and sum them (ln(L/); 
the alternative hypothesis). We also calculate the log likelihood 
obtained from the fitting for the whole period (ln(L const ); the 
null hypothesis). Then, TS var = -21 n(L const /L f ) is dist ributed 
as x 1 w i m three degrees of freedom (No lan et aDl2012l) . and 
we obtain ^ 2 /dof = 0.95. The probability for obtaining a value 
^f 2 /dof > 0.95 with three degrees of freedom is ~ 40%, which 
means that there is no indication for variability of Centaurus B 
at y-rays for the period spanned by the Fermi-L AT observations. 

3. Suzaku data 

We observed Centau rus B with the Suzaku X-ray satellite 
dMitsuda et al.l l2007h on 2011 July 16-18 (Observation ID 
806017010). The nominal pointing position was set to the core 
of the radio galaxy. Since the main goal of the analysis is to 
search for the diffuse X-ray emission from the radio lobes, here 
we fo cus on the data taken with the X-ray Imaging Spectrometer 
(XIS; lKoyama et al1l2007h . which consists of four CCD sensors 
each placed on the foc al plane of the X-ray Telescope (XRT: 
ISerlemitsos et al.ll2007l) . Currently two front-illuminated CCDs 
(XIS0 and XIS3) and one back-illuminated CCD (XIS1) are op- 
erational, and the observation was performed with all the XIS 
sensors set to full- window and normal clocking modes. In this 
paper we present the results obtained with XIS0 and XIS 3 (here- 
after FI). This is because the injection charge amount has been 
changed from 2 keV to 6 keV equivalent for the BI sensor on 
201 1 June 1, and the current non-X-ray background (NXB) esti- 
mation for XIS 1 using xisnxbgen does not reproduce accurately 
the actual level of NXB0. In the spectral analysis, we also ig- 
nore the data between 1.8 and 2.0 keV, because there exists a 
calibration uncertainty around the Si-K edgfl 

Since the radio lobes of Centaurus B are extended over al- 
most the whole Suzaku CCD chip, it is difficult to extract back- 
ground information from the data. Hence, we also analyze the 
Suz aku data for 4U 1344 -60 (Observation ID 705058010; see 
also ISvoboda et al.l (120121) ). which is located ~ 14' away from 
Centaurus B. We assume that the diffuse X-ray backgrounds are 
the same between these data given the proximity. This observa- 
tion was performed on 2011 January 11-14 with full- window 
and normal clocking modes. 

The screening criteria of XIS data are as follows. We uti- 
lized events with GRADE 0, 2, 3 ,4 and 6, then removed flicker- 
ing pixels by using cleansis. Good-time intervals were selected 
by the same criteria as described in The Suzaku Data Reduction 
Guid^ namely "SAA_HXD==0 && T_SAA_HXD>436 



&& ELV>5 && DYE_ELV>20 && ANG_DIST<1.5 && 
S0_DTRATE<3 && AOCU_HK_CNT3JNTML_P==l". In the 
screening and the following analysis, we used the HESoftB ver- 
sion 6.11 and calibration database (CALDBfl released on 2012 
February 11. The observation log and net exposures after the 
event screenings are summarized in Table 1 . 



3.1. Core 

Figure[5]shows the FI image in the 2- 10 keV band after subtract- 
ing the NXB and correcting for vignetting and exposure. We ex- 
tracted the events in a circular region with a radius of 4' centered 
at the Centaurus B core position and produced the XIS spectrum 
shown in Figure [6] The background region was selected from 
blank sky within the field of view (see Figure [5]). Response and 
ancillary response files were generated by xisrmfgen and xissi- 
marfgen, respectively. Then, we fitted the 0.6-9.0keV data with 
an absorbed power-law model (wabs*powerlaw) in XSPEC. The 
best-fit parameters are summarized in Table [2] 
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Fig. 5. Suzaku FI image (2 - lOkeV) of Centaurus B. The upper bright 
point-like source is the nucleus of the Centaurus B, and the lower source 
is 4U 1344-60. Green circle denotes the source extracted region for 
the Centaurus B core. Two green rectangles on the same CCD indicate 
the background regions. The diffuse X-ray background spectrum was 
produced from events within the lower white rectangular region. 



2 http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/node8.html 

3 http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/node8.html 

4 http://heasarc.nasa.gov/docs/suzaku/analysis/abc/ 



5 http : //heasarc . nasa. gov/lheasoft/ 

6 http://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/ 
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Table 2. Best-fit model parameters for the core of Centaurus B 
and the source N 



Parameter 
N H [xl0 22 



Core 



cm- 2 ] 



r x 

Normalization 



1.63 + 0.08 
1.62 + 0.05 



Source N 

1.55+0.56 

2.02+Hf 



Unabsorbed flux (5.1 + 0.4) x 10~ 12 
[erg cm -2 s _1 ] 



(l.l0!^)xl0-3 (L02^|) xlO- 4 



2.57 x 10" 



XiixVd.o.f) 



1.18 (366.13/310) 1.30 (23.37/18) 



Notes. (t) Af H and T x are absorption column density and photon index, 
respectively. The unabsorbed flux were integrated in 2-10 keV range. 
Errors are 90% confidence level. 



7 > 0.01 



.a u 
•a 




Energy (keV) 

Fig. 6. Suzaku/XIS spectrum of the Centaurus B core. XIS0 and XIS3 
spectra are shown in black and red, respectively. 



3.2. Diffuse Background 

A spectral analysis of faint diffuse emission requires the precise 
evaluation of the X-ray background emission. We extracted the 
background spectrum using the blank field of 4U 1344-60 (the 
IT x 2' rectangular region shown in Figure [5]). Here we con- 
sider contamination from 4U 1344-60 given a relatively large 
PSF (~ 2') of Suzaku XRT. FI count rate of the contamina- 
tion simulated by xissim ( Ishisaki et al was 0.0023 cts s 
in the 0.5 - 7.0 keV range, while the measured count rate 
was 0.01 14 cts s" 1 . We thus found that about 20% of the X- 
ray events within the considered rectangle are contaminated by 
4U 1344-60. Note that we have analyzed 14 ks Swift/XRT 
archival data and confirmed that there are no point sources in 
the rectangular region. We derive the 95% confidence level up- 
per limit for the unabsorbed flux of 1.1 x 10" 13 erg cm -2 s" 1 at 
0.5 - 10 keV, assuming a single power law with index 2.0 and 
N H = 1.05 x 10 22 cm" 2 . 

The diffuse X-ray background consists of three components: 
thermal emission from the Local Hot Bubble, a thermal com- 
ponent due to the Galactic halo, and a power-law component 
originating from a superposition of unresolved AGNs ("Cosmic 
X-ray Background", CXB). After subtracting both the NXB es- 
timated with xisnxbgen and the contaminant of 4U 1344-60, 
we fitted the background data in the 0.6 - 8.0 keV r ange with a 
model of mekal+wabs*(mekal+powerlaw) (see e.g. jYuasa et al.l 



Table 3. Best-fit parameters for the diffuse X-ray background 
spectrum derived from blank sky around 4U 1344-60 



Component 


Parameter 


Value 


T oral Hot Rnhhlp 

J— (VJCdl 11VJL DUUU1C 


kT TkeVl 


61 +0 09 

U ' Di -0.16 




z 


1.0 (fixed) 




Normalization 


(7.0 + 3.3) x 10- 4 




N H [xl0 22 cm" 2 ] 


1.06 (fixed) 


Galactic halo 


kT (keV) 


17 +009 

U>1 7 -0.06 




Z 


1.0 (fixed) 




Normalization 


26 +338 

^^"-0.24 


CXB 


r x 

Normalization 


9 99+0.23 
^'^ -0.22 

|) x 10 - 3 




xWld-o.f) 


1.02 (63.41/62) 



Notes. (t) kT and Z are temperature and metal abundance in terms of the 
solar value of the plasma, respectively. Errors are 90% confidence level. 



120091) . In the fit, the absorption column density was fixed to the 
value in that direction^ (N n = 1.06 x 10 22 cm" 2 ). The value of 
A/r obtained in the spectral analysis of the core is larger than this 
value possibly due to the absorption of the host galaxy. The best- 
fit spectra and parameters are shown in Figure [7] and Table [3j 
respectively. The derived power-law index does not agree with 
the "standard" CXB index of 1.41, but is sim ilar to that of the 
Galact ic ridge emission (index -2.1, see e.g. jRevnivtsev et all 
(120061) ). This could be explained by the fact that the analyzed 
field is located close to the Galactic plane. 



o.oi 




Energy (keV) 

Fig. 7. Suzaku/XIS FI spectrum of the X-ray background obtained from 
blank sky near 4U 1344-60. Dashed, dot-dashed, and dotted lines rep- 
resent the components of the Local Hot Bubble, the Galactic halo, and 
the CXB, respectively. 



3.3. Lobes 

We found two faint point sources within the radio contours of the 
lobes of Centaurus B in the 0.5 - 2 keV and 2-10 keV Suzaku 
FI NXB -subtracted exposure-corrected images (see Figure [8]). 
After masking the core and these two X-ray sources with 4' and 



http://heasarc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 
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Fig. 8. (Left) Suzaku FI 0.5 - 2keV image of Centaurus B. (Right) 
Suzaku FI 2 - lOkeV image of Centaurus B. The central bright source 
is the nucleus of the sy stem. Green conto urs denote the 843 MHz radio 
intensity obtained from McAdam| dl991h . Magenta ellipse in the image 
denotes the source extracted region for the lobes. The three magenta 
circular regions containing X-ray point sources are excluded from the 
analysis of the emission from the lobes. 



2' radii circles, the elliptical region (22' x 10' in size) was ex- 
tracted for the spectral analysis of the lobes. NXB was estimated 
with xisnxbgen. The contamination from the Centaurus B core 
was simulated by xissim, where the values in Table [2] were used 
for the spectra of the core. The spectrum of diffuse X-ray back- 
ground of this region was simulated with XSPEC by using the 
values tabulated in Table [3] Figure [9] shows the observed spec- 
trum together with the three background components. We uti- 
lized xissimarfgen to produce the ancillary response file assum- 
ing a uniform surface brightness in the elliptical region. 




Energy (keV) 

Fig. 9. Suzaku/XIS FI spectrum of the Centaurus B lobe region. The ob- 
served spectrum are shown in black. Green, red and blue spectra denote 
the background components: the diffuse X-ray background, the NXB, 
and the contaminant from the core, respectively. 



As shown in Table 3, the ~25% statistical error in normalization 
dominates the ~3% systema tic error in the estimation of NXB 
spectrum (Tawa et al. 2008) up to about 6 keV (see Figure [7]). 
Hence, we increased/decreased the normalization of power-law 
component by 25% and produced faked X-ray background spec- 
trum. Then after subtracting it we searched for the residual emis- 
sion above 2 keV. When we increased the normalization by 
25%, the residual spectrum does not show significant power-law 
component and the best-fit normalization corresponded to zero 
(rx = 1 .7 is assumed following the radio data; see §@]). To derive 
the upper limit for the lobe X-ray emission, we considered the 
most conservative case, namely the normalization of the power- 
law component of X-ray background decreased by 25%. Then 
we fitted the residual spectrum by assuming an absorbed power- 
law shape with only normalization set free and the other parame- 
ters fixed as r x = 1.7 and N n = 1.06 x 10 22 cm" 2 . 9 The obtained 
90% confidence level upper limit for the unabsorbed 2-10 keV 
emission of the Centaurus B lobes is < 8.7 x 10" 13 erg cm -2 s" 1 . 

As shown in Figure [8] one of the masked sources (source N), 
located at (a, <5)=(20ffi721, -60?506), has a hard X-ray spec- 
trum. This source is positionally coincident with the diffuse-like 
emission seen in the ASCA GIS 1.5 - 3.0 keV image within th e 
South- Western radio lobe (see Figure 3 in lTashiro et al.|[l998h . 
To derive spectral parameters of the source N, we extracted a 
circular region with a radius of V centered at this source, and 
produced a FI spectrum as shown in Figure [TO] The background 
region was selected from the blank field of the same CCD chip. 
The spectrum was well represented by the absorbed power-law 
model as shown in Table [2] We also searched for possible coun- 
terparts of source N in other wavelength using SIMBAE0, and 
found an infrared source just at the position of source N. This co- 
incidence again supports the possibility that the point-like X-ray 
emission is unrelated to the radio lobe. 



io-- 



io- 4 




Energy (keV) 

Fig. 10. Suzaku/XIS FI spectrum of the source N. 



Since the X-ray events of the diffuse background emission 
within the selected background region are small, statistical errors 
of the derived parameters are relatively large (see Table [3]). In 
particular, the error in the normalization of the power- law com- 
ponent of the X-ray background significantly affects the evalu- 
ation of the faint diffuse X-ray emission from the radio lobes. 



4. Radio data 

An 843 MHz image (iMcAdaml Il99lh from the Molonglo 
Observatory Synthesis Telescope (MOST), with 43 " resolution 
and excellent surface brightness sensitivity, provides a reliable 
total integrated flux density of ~ 150 Jy. The northeastern lobe is 

8 http ://simbad. u- strasbg . fr/simbad/ 



7 



J. Katsuta et al.: Fermi-L AT and Suzaku Observations of the Radio Galaxy Centaurus B 



substantially brighter with an integr ated flux of - 93 Jy while the 
southwestern lobe contains ~ 57 Jy. Jone s et al ] (120011) used inte- 
grated flux densities from the literature, ranging from 30 MHz to 
5 GHz, to derive a radio spectral index of a = 0.73. As noted in 
§[T] the 1 .4 GHz radio power puts Centaurus B near the boundary 
of the FRI/FRII luminosity division. The bright jet and diffuse 
lobe emission resembles the morphology of an FR I radio source, 
while the northeastern lobe has a shar p edge reminiscent of an 
FR II radio galaxy (see in this context iGopal-Kris hna & Wiital 
l200ttlLaingetalJl2011h . 

Jones et al. (120011) also present observations of the entire 
source using the Australian Telescope Compact Array (ATCA) 
at 4.8 and 8.6 GHz. We are unable to obtain reliable integrated 
flux densities for the lobes from these images because of the 
large angular extent of the source O 24'). A substantial amount 
of flux is lost, since the shortest baseline between antennas 
(~ 20 m for mosaicked images from the ATCA) is too large 
to recover the emission on the largest spa tial scales. The entire 
sourc e has also been imaged at 20 GHz (iBurke-Spolaor et all 
2009]). A gain, the integrated fluxes cannot be used in our SED 
modeling due to the missing short spacings, although the results 
provide an accurate representation of the core, jet and hotspot 
structures. Also, in this particular case, potential imaging arti- 
facts due to the bright core may contaminate the faint lobe emis- 
sion at 20 GHz. 

Since our line of sight to the source passes within 1.7° of 
the Galactic plane, single dish measurements can be contam- 
inated by the Galactic emission (particularly at low radio fre- 
quencies). The Parke s-MIT-NRAO survey conducted at 4.8 GHz 
(iCondon et al.lll993l) provides us with a second suitable image 
for extraction of integrated fluxes. The Parkes 64 m dish im- 
age has high enough angular resolution (~ 5') to resolve the 
source into two lobes, with integrated flux densities of ^ 29 Jy 
and 18Jy for the northeastern and southwestern lobes, re- 
spectively. Furthermore, the observing frequency is also high 
enough, at this spatial resolution, to avoid any substantial con- 
tamination from diffuse Galactic emission. The core an d inner jet 
flux contributes ~3 Jy at 5 GHz dMarshall et al.ll2005l) . which is 
a small fraction of the total integrated emission of ~47 Jy. Since 
the lobes are brighter at lower frequencies, the contribution from 
the core becomes negligible at frequencies below 5 GHz. Hence, 
the integrated fluxes from 30 MHz to 5 GHz provide accurate 
representations of the integrated emission from the lobes. 

The newly launched Planck satellite ( Planck Collaboration 
et al. l201 lal) has observed the sky in nine frequency bands from 
30-857 GHz with the first all- sky c ompact- source catalogue of 
high-re liability sources presented in iPlanck Collaboration et al.l 
(1201 lbl) . They used detailed source extraction algorithms to re- 
liably detect discrete sources immersed in a diffuse foreground 
(e.g., the Galactic plane). From their catalogue, we were able 
to extract integrated fluxes for Centaurus B at 30, 44, 70, 100 
and 143 GHz. The spatial resolution of 33' at 30 GHz means 
that Centaurus B is unresolved and the total flux density accu- 
rately represents the entire flux (lobes plus core) from the source 
at this frequency. We cannot discriminate between the contri- 
bution of the core and lobe fluxes from the available data at 
higher radio frequencies, however, we expect the flat- spectrum 
core to become increasingly dominant over the steep- spectrum 
lobes above 30 GHz. The extracted 44 GHz flux density appears 
to slightly underestimate the total flux of the source at this fre- 
quency (based on power-law fit to all Planck data points), even 
though the spatial resolution of 27 ' should be sufficient to in- 
clude the entire emission from the source. The spatial resolution 
of Planck at 70, 100 and 143 GHz of 13', 10' and T means that a 



large fraction of the flux from the lobes of Centaurus B is miss- 
ing. 

5. Discussion 

The observed GeV emission of Centaurus B could be produced 
within the unresolved core of the system, in the analogy to 
some other radio galaxies detected by the Fermi-LAI such as 
NGC 1275 jXbdo et al.l 12009^ Kataoka_eL_al] HoH, M87 



dAbdo et al.ll2009bl:lAbramowski et al.ll2012l) . 3C 120 or 3C 1 1 1 
(iKataoka et al.ll201ll) . In those cases, the broad-band nuclear 
spectra as well as the established y-ray variabilities indicate 
that the y-ray production is dominated by the inner (pc/sub-pc 
scale) relativistic jets observed at intermediate viewing angles 
("misaligned blazars"). The jet/counterjet radio brightness and 
polarization asymmetries on slightly larger scales in the source 
studied here imply intermediate jet inclinations and relativistic 
beaming involved, supporting the idea that the nucleus may con- 
tribute to the observed y-ray emission, at least at some level. 
On the other hand, the LAT analysis of Centaurus B provides a 
weak hint of the spatial extension of the y-ray counterpart, as 
well as no indication of significant flux variability in the GeV 
range, which suggest that the extended (100-kpc scale) lobes in 
the system could be considered as a promising emission site of 
the high-energy photons as well. The arguments based on the 
morphology of the analyzed y-ray emitter and its timing proper- 
ties are nonetheless rather weak, due to the limited photon statis- 
tics. 

At X-ray frequencies, the core of Centaurus B appears as 
a relatively b r ight s ource. In their analysis of the ASCA data, 
iTashiro et al.l (1 19981) found the 2 - lOkeV nuclear unabsorbed 
flux of ^ (6.4 ± 0.2) x 10" 12 erg cm -2 s" 1 with a power-law model 
ofT x - 1.64 ±0.07 and N n ^ (1.76 ± 0.1 1) x 10 22 cm" 2 possibly 
in excess over the Galactic val ue. This is in agreeme nt with the 
Chandra snapshot obtained by Marsha flet al. I (12005). who mea- 
sured the 1 keV nuclear flux 870+ 24 q nJy and the best-fit power- 
law model with Tx - 1 .63 ±0. 17. Our new Suzaku exposure con- 
firms the previous studies, revealing the unabsorbed 2 - lOkeV 
flux of the core ^ (5.1 ± 0.4) x 10" 12 erg cm -2 s" 1 marginally 
lower than that implied by the ASCA observations but con sistent 
with that detected by XMM-Newton (ITashiro et al.ll2005l) . with 
T x - 1.62 ± 0.05 and N n - (1.63 ± 0.08) x 10 22 cm" 2 . The un- 
absorbed X-ray (2 - lOkeV) luminosity of the core, Lx - 2 x 
10 42 erg s" 1 , is roughly comparable to the y-ray (0.2 - 200 GeV) 
luminosity of the Centaurus B system, L y 6 x 10 42 erg s" 1 (see 
Figure[TT]). m addition, the inner parts of the main (SW) jet were 
resolved at arcsec scale with Chandra dMarshall et al.l2005l) . The 
detection of the SW jet on larger scales was also claimed in the 
preliminary analysis of the XMM data (iTashiro et al .1120051) . If 
the observed X-ray continuum of the nucleus corresponds to the 
beamed non-thermal (inverse-Compton) emission of the unre- 
solved jet, the observed y-ray flux could be explained well as a 
high-energy tail of the same nuclear emission component. 

A caveat is that in radio galaxies the observed X-ray nu- 
clear emission is typically dominated by the accretion disks and 
disk coronae, and not by the misaligned unresolved jets. This 
seems to be particularly robust in the case of powerful objects 
accreting at high rates, like Broad-Line Radio Galaxies fo r ex- 
ample (e.g., lEracleous et al.l l2000l : IKataoka et al.ll201ll) . For 
the low-power objects accreting at low rates, such as FRI ra- 
dio galaxies, their entire observed X-ray continua may be due 
to th e innermost parts of the jets (iDonato et al.ll2004l; Evans 
et al. 120061 : lHardcastle et al.ll2009bl) . This controversy relates 
in particular to low-power radio galaxies detected by the Fermi- 
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LAT (see the discussion in Abdo et al ]l2009a.b. 2010a, and ref- 
erences therein). Centaurus B seems to be an intermediate case 
between low- and high-power sources, and hence one may ex- 
pect a mixture of the nuclear jet and accretion-related emission 
components at X-ray frequencies; the observed spectral proper- 
ties of the X-ray core in the system suggest however a dominant 
contribution from the jet. 

In order to test the hypothesis that the observed y-ray emis- 
sion of Centaurus B originates in the innermost parts of a mis- 
aligned relativistic jet, we apply the synchrotron self-Compton 
(SSC) model to the collected broad-band data of the unresolved 
core in the systenfl assuming that the emission site can be ap- 
proximated as a homogeneous spherical region with radius R\ 
filled uniformly with tangled magnetic field and radiating elec- 
trons. This is the first time that the SSC model is applied to the 
misaligned relativistic jet in Centaurus B and compared with the 
compiled multi wavelength dataset for the unresolved core. As 
shown in Figure [TT] the model successfully explains the data 
for the following parameters: jet radius R' = 0.03 pc, jet mag- 
netic field B = 0.095 G, and jet Doppler factor 6 = 2. For 
the radiating electrons we assume a power-law energy spectrum 
dN' e (y)/dy oc y~ 2 - 2 , where y is the Lorentz factor, between ener- 
gies ymin = 30 and y^ v = 2 x 10 3 , breaking to dN' e (y)/dy oc y~ 32 
above y^ Y and continuing as such up to y max = 10 5 . The fit re- 
quires normalization of the energy spectrum such that the ra- 
tio of the comoving electron and magnetic energy densities is 
^eq = U'elU'B = AH these parameters are broadly consistent 
with the ones inferred for the cor es of othe r GeV-detected radio 
galaxies (see lAbdo et al. 2009aHbl. feOlOal) . Also the jet power 
emerging from the model presented here should be considered 
as reasonable: assuming no significant jet radial structure, the 
model value of the jet Doppler factor 5 = 2 implies jet bulk 
Lorentz factors Tj ■ = 4 - 7 for the jet viewing angles within the 
range 0j = 20° - 25°, and hence the jet luminosity, roughly, Lj 
nR f2 cT 2 U' e ~ 10 45 ergs _1 . This estimate neglects likely contri- 
bution of protons to the jet kinetic flux, and therefore should 
be considered as a lower limit only. The implied jet-counterjet 
brightness asymmetry, R = [(1 + /3jCOs6j)/(l -J3j cos 6 j)] £ = 

(2TjS - 1) £ , where J3j = ^jl - Tt 2 is the jet bulk velocity and 

s = 2 + a assuming a steady jet and the emission spectral in- 
dex a, is on the other hand rather high, R > 10 3 f or a > 0.5, 
much greater than the one measured on larger scales (iJones et al.l 
l200lh . This is however consistent with the established idea of 
the jet deceleration from sub-pc to kpc scales — the process in- 
evitably accompanied by the formation of radial stratification of 
the jet — in low and intermediate power radio galaxies (FRIs 
and FRI/FR lis) such as Centaur us B ( see in this context, e.g., 
lLaing et al.|[l999l : lLaing & Bridlell2002h . 

Let us now focus on the extended lobes in the Centaurus B 
system. For modeling the broad-band SED of the lobes we use 
the the t otal radio fluxes i n the range 30 MHz -5 GHz as pro- 
vided by lJones et al.l (l200lh . as well as the 30 GHz flux measured 
by Planck, which should represent the integrated emission of the 
lobes (see §|4]). We also take into account the X-ray upper limit 
of our Suzaku observations, ^-lokev < 8.7 x 10" 13 erg cm -2 s" 1 
assuming T x = 1.7 as shown in § [3] We note that this upper 
limit is only marginally consistent with the flux found by Tashiro 




Fig. 11. Broad-band SEDs of Centaurus B and the model curves for the 
emission of the extended lobes and core in the system. Re d circles de- 
note t he total radio fluxes within the range 30 MHz -5 GHz (IJones et aTl 
1200 lL and references therein) and at 30 GHz measured by Planck. The 
radio fluxes for the unresolved core, p lotte d as blue stars, are taken from 
IJones et all fcOOlh . lFev et all d2004l) and lBurke-Spolaor et all (120091) . 
The ASCA measurements of t he lobes' and core spe ctra within the range 
0.7 - lOkeV, as discussed in lTashiro et al.1 (119981) . are denoted by red 
and blue bow-ties, respectively. The Chandra 0.5 - 7 k eV spectrum of 
the un resolved core is given by cyan bow-tie following Ma rshall et al.1 
d2005h . The Fermi-LAI fluxes represented by magenta squares corre- 
spond to the analysis presented in §[2]of this paper. The black bow-tie 
denotes the 2 - 10 keV emission of the nucleus emerging from the anal- 
ysis of our Suzaku data, while the thick black line shows the 2 - 10 keV 
Suzaku upper limit foe the lobes' emission (see § [3j. Thick gray curves 
denote the broad-band emission of the lobes modeled in in the frame- 
work of the synchrotron and inverse-Compton scenario, while black 
dashed curve corresponds to the SSC model of the unresolved core in 
the s ystem, as discussed in §(5] 



9 Note that as the observed core radio fluxes correspond to a super- 
position of different self-absorbed jet components, in the fitting proce- 
dure we consider these as upper limits for the emission of the modeled 
emission site, for which the rest-frame turnover frequency related to the 
synchrotron self-absorption (SSA) process reads as vj sa = 9 GHz. 



et al. (Il998h . Fi.mv^r ^ (9.8 ± 3.8) x 10" 13 erg cm" 2 s" 1 . We 
assume that the radio emission is due to the synchrotron radi- 
ation of non-thermal electrons distributed isotropically within 
the homogeneous lobes. We approximate the lobes for sim- 
plicity as two spheres with radii R ^ 100 kpc, so that the to- 
tal volume of the emission region is Vi ^ 2 X (4n/3)R 3 . The 
high-energy continuum, on the other hand, including the X- 
ray and y-ray fluxes, is ascribed to the inverse-Compton (IC) 
emission of the same electron population. As in the case of 
Centaurus A (I Abdo et al.ll2010bl) . we consider three relevant tar- 
get photon fields: the CMB, the extragalactic background light 
(EBL) photons within the infrared-to-optical range, and the 
starlight of the host galaxy In ou r calculations we use the EBL 
model of iRaue & Mazinl (120081) . and note that this particular 
choice does not affect the obtained results. For the starlight of 
the host galaxy, we use a template for the giant elliptical nor- 
malized toJheAA_><l&^ flux of 
0.42 Jy (lLaustsen et al.l 1 19771 : [Schroder et al.ll2007l) . This gives 
a V-band luminosity for the host of Ly - 10 44 ergs -1 , and the 
corresponding energy density of the starlight within the lobes 
^star - L star /4nR 2 c ~ 3L Y /4nR 2 c < 10 -14 erg cm" 3 . The results 
of the modeling are presented in Figure [TT] 

As shown in Figure [TT] the simplified model can explain 
the radio and y-ray fluxes well, not violating the X-ray upper 
limit. For the model parameters, the IC spectrum is dominated 
by the CMB component, with negligible contribution from the 
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EBL and the starlight of the host galaxy. A broken power-law 
electron energy distribution dN e /dy is set to oc y~ 2A for y < ybr 
and dN e /dy oc y~ 3A for y > y^ Y , with minimum, break, and 
maximum Lorentz factors of y^n - 1, Tbr = 3.5 x 10 5 , and 
Tmax = 1.5 x 10 6 , respectively. The lobes appear to be slightly 
out of the energy equipartition, with the corresponding electron- 
to-magnetic field energy density ratio 7/ eq = U e /UB ^4, and the 
mag netic field intensity B ^ 3.5 uG (see in this context T ashiro 
et al.Tl998; Katao ka & Stawarzll2005l: ICroston et al.ll2QQ5l: Isobe 
et al. l201ll) . The lobe pressure stored in ultrarelativistic elec- 
trons and magnetic field is /? e +B - 8 x 10" 13 dyncm -2 , while the 
total energy injected by the jets over the source lifetime (a sum 
of the work done in displacing a volume Vt of surrounding gas 
at pressure /? e +B and the energy of the material inside the cav- 
ity) is Etot - 4 p Q+ B Vi ~ 8 x 10 59 erg. These values, derived 
from the crude approximations in the mod eling, are typical for 
intermediate-power radio galaxies (s ee e.g., ICroston et al.ll2005l ; 
iBirzan eta! .2008; Isobe et al. 2009). 

It can be seen from Figure [TT] that the model fit depends cru- 
cially on the high-energy tail of the electron energy distribution 
within the lobes, which is currently unconstrained by the data. 
To account for the Fermi-LAI fluxes, the electron break Lorentz 
factors within the lobes must be as high as y^ = 3.5 x 10 5 , giving 
the expected break synchrotron frequency Vb r - 1 .8 x 10 12 Hz for 
B ^ 3.5 //G. This would then correspond to the spectral age of 
the lobes 



T rad - 50.3 



(5/10//G) 



1/2 



(£/10yUG) 2 + 0.1 



I Vbr r 1/2 ^ 
(GHz") My " 



3Myr, 



(1) 



according to the standard synchrotron a ging theory (see, e.g., 
iKonar et al.ll2006l : lMachalski et al.ll2007l) . This is a rather young 
age for the observed ~ 400 kpc extension of the entire radio 
structure, implying either particularly high advance velocity of 
the jets in the system, v a dv <; 0.1 c in order for the dynamical 
age Tdyn ~ 27?/ Vadv to be comparable with r ra d, or efficient re- 
acceleration of ultrarelativistic particles stored within the lobes. 
Such a re-acceleration process (the exact nature of which we do 
not specify) would result in a younger spectral appearance of the 
lobes. In this case, the age derived based on the spectral anal- 
ysis, T ra d ~ 3 Myr, could be much shorter than the true age of 
the system, namely Td yn ~ 2R/v a ^ > 10 Myr with the jet ad- 
vance velocities v a dv < 0.1 c, which are typically derived for the 
evolved radio galaxies (linear sizes > 10 kpc; see the discus- 
sions in Scheuer 1995; Alexander 2000; Machalski et al. 2007; 
iKawakatu et al.ll2008l) . 

Here we favor the re- acceleration scenario, since the alter- 
native explanation involves uncomfortably high values of v a dv ~ 
0. 1 c, exceeding by orders of magnitude the sound velocity in the 
surrounding medium c s = ^j5kT/3m p ~ 0.001 (r/10 7 K) 1/2 c 
(with the expected T ~ 10 7 K). This would imply the pres- 
ence of a strong bow- shock driven in the ambient medium 
by the supersonically expanding lobes, with the Mach number 
M = Vadv/Cs » 10 unprecedentedly high as for an evolved ra- 
dio galaxy; we emphasize that all the bow- shocks discovered till 
now in such systems are c haracterized by very low M ach num- 
bers, typically At ~ 2 (see lSiemiginowska et all2012l and refer- 
ences therein), with the prominent exception of the At - 8 shock 
detected around the Southern inner lobe in the Centaurus A sys - 
tem, albeit on much smaller (kpc) scales (ICroston et al.ll2009h . 
We also note that similar conclusions regarding particle re- 
acceleration in the extended lobes have been presented in the 
case of two other radio galaxi es resolved by t he Fermi-LAI, 
Centaurus A (lAbdo et alJ2010bh and NGC 6251 (Takeuchi et al. 



120121) . Clearly, high-frequency (> 30 GHz) radio observations 
by ground-based interferometers with the appropriate resolution 
and sufficiently short baselines are required to make further ro- 
bust statements in this respect, since the association of the y-ray 
emission site in Centaurus B could be confirmed further by con- 
straining the high-energy tail of the lobes' synchrotron contin- 
uum. 

It is also important to determine the level of the diffuse non- 
thermal X-ray emission in Centaurus B by ultra-deep exposures 
with modern X-ray telescopes. If the diffuse non-thermal X-ray 
emission is significantly below the upper limits derived in this 
paper, then the association of the y-ray emitter with the radio 
lobes should be viewed as unlikely. Such a low level of the X- 
ray non-thermal halo would imply that the lobes are dominated 
by magnetic pressure, contrary to previous claims presented in 
the literature, or at least 7/ eq - 1 . This is because the X-ray-to- 
radio luminosity ratio is determined by the magnetic field energy 
density (or pressure) within the lobes, Lx/Lr = U cmrj /UB oc B~ 2 , 
where t/ cm b is the energy density of the CMB photon field. The 
lower X-ray flux for a given radio flux implies a lower value of 
the 7/eq parameter. In particular, noting that the inverse-Compton 
X-ray luminosity Lx oc £/ cmb x U e oc 7/ eq £/b, one may find 
7/eq oc L 2 for the fixed Lr and U cm \y. In the modeling presented 
above, the expected X-ray luminosity of the lobes is only a factor 
of ^ 1.4 below the provided Suzaku upper limit, and the corre- 
sponding equipartition parameter is 7/ eq = 4. On the other hand, 
if the X-ray flux is a factor of 3 below the current upper limit, the 
equipartition parameter is rj Qq < 1 and magnetic pressure dom- 
inates over the pressure of ultrarelativistic electrons within the 
lobes of CentaurusB. 



6. Conclusions 

In this paper, we present a detailed analysis of about 43 months 
accumulation of Fermi-LAI data and of newly acquired Suzaku 
X-ray data for CentaurusB. We detect the source in GeV pho- 
ton energies at 4.5 <x significance (TS = 29). The implied y-ray 
(0.2 - 200 GeV) luminosity of the system is L y 6 x 10 42 erg s" 1 , 
and the observed continuum within the LAT range is best mod- 
eled by a power-law function with the photon index Y y = 2.6 ± 
0.2. At X-ray frequencies we detect the core of Centaurus B with 
Suzaku at the flux level consistent with the previous measure- 
ments with ASCA, XMM-Newton and Chandra. The unabsorbed 
2 - lOkeV luminosity of the core is Lx - 2 x 10 42 erg s" 1 , and 
the spectrum is best fitted with the power-law model with pho- 
ton index Tx - 1.62 ± 0.05. We do not detect any significant dif- 
fuse X-ray emission of the lobes, with the provided upper limit 
^2-iokeV < 8.7 x 10" 13 erg cm -2 s" 1 , assuming the photon index 
Tx = 1.7 . The upper limit is marginally consistent with the flux 
found by lTashiro et al.l (E998) in their analysis of the ASCA data, 
^2-iokeV - (9.8 + 3.8) x 10" 13 erg cm" 2 s _1 . 

By means of broad-band modeling we show that the ob- 
served y-ray flux of the source may in principle be produced 
within the lobes via the inverse-Compton scattering of the CMB 
photon field by radio-emitting electrons, if the diffuse non- 
thermal X-ray emission component is not significantly below 
the derived Suzaku upper limit. This association would imply 
efficient in-situ acceleration of ultrarelativistic particles within 
the lobes dominated by the particle pressure, in direct anal- 
ogy to two ot her radio galaxies resolved by Fer mi-LAI, namely 
Centa urus A (lAbdo et al.ll2010bl) and NGC 625 1 dTakeuchi et al J 
120121) . However, if the diffuse X-ray emission in the studied ob- 
ject is much below the derived Suzaku upper limits, then the ob- 
served y-ray flux is not likely to be produced within the lobes, 
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but instead within the unresolved core of the Centaurus B ra- 
dio galaxy. Assuming the one-zone synchrotron self-Compton 
model, we show that this possibility is also justified by the broad- 
band data collected for the unresolved core of Centaurus B, in- 
cluding the newly derived Suzaku spectrum, and the expected pa- 
rameters of the misaligned nuclear relativistic jet. Interestingly, 
in such a case the extended lobes in the system could be domi- 
nated by the pressure of the magnetic field. 
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